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Abstract. Using the theory of optimal rocket trajectories in general relativ- 
ity, recently developed in IHNlll . we present a candidate for the minimum 
total integrated acceleration closed timelike curve in the Godel universe, and 
give evidence for its minimality. The total integrated acceleration of this curve 
is lower than Malament's conjectured value |Mal84] . as was already implicit in 
the work of Manchak IManlll : however, Malament's conjecture does seem to 
hold for periodic closed timelike curves. 



Introduction 

It is well known that the Godel universe contains closed timelike curves |God49al 
IHE95] . These curves are not geodesies, and so any particle following them must 
accelerate. Godel himself remarked that if the acceleration is provided by a rocket 
then the required amount of fuel is prohibitively high |God49bl IMal871 IOS03] . This 
amount can be computed from the well known rocket equation |Ack46[ IHNll) . 
which relates the initial and final rest masses of the rocket: 



miji) = TO(ro)exp / a(T)dT 
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where r is the proper time, m(T) is the rest mass, a(r) is the magnitude of the 
proper acceleration and v is the ejection speed. Since v must be smaller than the 
speed of lighlQ, i; < 1, we have 

TO(ro) 



TO(ri) 

where 



> exp(TA), 



TA= f a{T)dT 

J To 



is the so-called total integrated acceleration. Malament |Mal841 IMal85| proved 
that TA > ln(2 + %/5) ~ 1.4436 for any closed timelike curve in the Godel universe, 
and conjectured that in fact TA > 9 + 6\/3 ~ 27.6691. Notice that if this 
conjecture is true then 

> 10^ 



m(ri) 

for any rocket following a closed timelike curve in the Godel universe! Manchak 
[Manllj . however, proved that there exist closed timelike curves in the Godel uni- 
verse with TA < 27rv/9 + 6^3. These curves are not periodic, in the sense that the 
initial four-velocity is not equal to the final four-velocity, and it seems likely that if 
one includes the integrated acceleration necessary to boost the final four- velocity to 



Partially supported by FCT (Portugal). 

"'^We use the conventions of |MTW73] , including geometrized units (c = G = 1). 
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the initial four-velocity then Malamcnt's conjecture still holds. Nevertheless, Man- 
chak did succeed in showing that one can have time travel in the Godel universe 
with total integrated acceleration smaller than Malament's bound. 

The situation is analogous to the problem of finding a closed curve in Euclidean 
space for which the total integrated geodesic curvature with respect to the arclength 
(the Riemannian analogue of the Lorentzian total integrated acceleration) is min- 
imum. If we require the curve to be periodic, then the minimum is 2tt, and it is 
attained by a circle, for example. If, however, we do not require the initial and final 
directions of the curve to match up, there is no minimum; the infimum is tt, and is 
approached by a sequence of closed curves consisting of an arc of circle subtending 
an angle tt -f i (for n £ N) and two line segments tangent to its endpoints, which 
therefore meet at an angle i (see Figured]). The intersection of the line segments is 
meant to be the curve's initial and final point. Notice that the integrated geodesic 
curvature required for "turning the corner" at this point (thus making the curve 
periodic) is precisely tt — i. 



Figure 1 . Sequence of closed curves with total integrated geodesic 
curvature approaching tt. 

In this paper we present a candidate for the minimum total integrated accelera- 
tion closed timelike curve in the Godel universe, and give evidence for its minimal- 
ity. The total integrated acceleration of this curve turns out to be approximately 
24.9927. If one includes the integrated acceleration necessary to boost the final four- 
velocity to the initial four- velocity, however, total integrated acceleration increases 
to about 28.6085, and so Malament's conjecture is not violated. 

The candidate curve is more complicated than one would guess from the Eu- 
clidean space analogy, and is obtained by applying the general theory of optimal 
rocket trajectories in general relativity developed in f HNll) . For the reader's con- 
venience, this theory is summarized in section[l] and the main features of the Godel 
universe are reviewed in section [2j In section |3] we solve the Jacobi equation for 
a timelike geodesic, and in section 2] we compute the circular optimal accelerated 
arcs, both of which are used in section [5] to construct the candidate curve and give 
evidence for its minimality. 



In [HNllj the classical Newtonian theory of optimal rocket trajectories devel- 
oped in |Law63j was generalized to the general relativity setting. There it is shown 




1. The rocket problem in general relativity 
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that optimal trajectories are expected to be continuous, sectionally smooth time- 
like curves, obtained by piecing together free-fall (geodesic) and accelerated arcs, 
possibly with instantaneous (Dirac delta) accelerations at the junction points. If 
the trajectory minimizes the total integrated acceleration then it satisfies the dif- 
ferential equations 

(1) < V[/P^ = -q^ + aJ7^ 

where U is the four-velocity, a is the magnitude of the proper acceleration, P is 
an auxiliary vector field, called the primer, and Rafi^iu are the components of the 
Riemann curvature tensor. All quantities above are continuous at a junction point 
except possibly U and P, which change by a boost with positive parametei0 in 
their common plane if there is an instantaneous acceleration. The magnitude of p 
the primer, however, is at the junction points, and satisfies p < 1 on free-fall 
arcs and p = 1 on accelerated arcs and instantaneous accelerations. Moreover, the 
equations above admit the first integrals 

(2) P^C/'^ = q^,U^' = 0. 

Finally, it is easily seen from |HN11) that if the initial or final four-velocities are 
not specified then we must have P = at those points. 

For free-fall arcs (a = 0), conditions ([T]) reduce to the geodesic and Jacobi 
(geodesic deviation) equations, with the first integrals ([2]) requiring the Jacobi field 
P to be orthogonal to the four-velocity U . 

For accelerated arcs, we have from and ([2]) 

q^C/" = ^ (V[/(Z^)!7'' + q^.{VuU'') = ^ g^P^ = 0. 
Differentiating this new first integral yields 

(3) (Vc/gjP^ + <z^(Vc/P'') = ^ q^g" = 
where 

K = -Ra,f}^,P''U^P''U'' 

is the sectional curvature of the tangent plane spanned by U and P. The magnitude 
a(T) of the proper acceleration must then be chosen so that ^ holds at each instant. 
Using ([T]) and ©, we can easily rewrite ^ as 

(4) (Vc/P^)(Vc;P^)+a2 = -if. 

2. The Godel universe 

The Godel universe is a space-time homogeneous solution of the Einstein equa- 
tions with a pressureless perfect fluid and a negative cosmological constant. It 
corresponds to the Lorentzian metric in given in cylindrical coordinates by 



ds^ = - (dt + V2{1 - coshr)dip] + dr^ + sinh^ r dip'^ + d 



^This boost parameter is exactly the contribution of the instantaneous acceleration to the total 
integrated acceleration 
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where we have chosen unit^ such that the cosmological constant is A = — ^ and 
the matter density is p = (see for instanc^^ |OS03| ). The fluid worldlines are 
tangent to the timelike Killing vector field ^ , and have vorticity given in our units 

by a; = (so the fluid elements are rotating with angular velocity about 

the z-axis with respect to the local compass of inertia). 

In what follows we will restrict ourselves to a totally geodesic submanifold of 
constant z (which we will still call the Godel universe), with metric 

ds^ =^ - (dt + \/2{l - cosh r)d'f) + dr^ + sinh^ r dip^ 



(the nontrivial part of the Godel metric) . Notice that the quotient of this manifold 

d_ 

dt 



by the isometric K-action generated by ^ is the hyperbolic plane of radius 1 (hence 



our choice of units). 

The Godel universe famously contains closed timelike curves (for instance the 
integral curves of ^ for r sufficiently large). For any timelike curve with proper 
time T we have ^ 

-dr'^ = -{dt + \/2{l-coshr)dipj + df , 

where dt^ is the line element of the hyperbolic plane. Solving for dt yields 



dt = \/2(coshr - l)dLp + ^/dP + dP^ 

for future-directed curves (corresponding to the positive square root). Therefore if 
the projection 7 of the timelike curve on the hyperbolic plane is closed then the 
coordinate time difference between its endpoints is 

At=(hdt = V2(l (cosh 7 

J J y J ^ 

It is easy to check that the velocity v of the timelike curve with respect to the 
matter satisfies 

dl V /, /dr \ ^ 1 




dr y'l — 1;2 y \dl J V 
If the curve 7 is simple and bounds a region D then we can use the Stokes theorem 
to write 

(5) At ^ V2 [ [ sinhrdr Adip + (f — = V2A+ (f —, 

JJd J-y V JjV 

where A is the oriented hyperbolic area of D (negative or positive according to 
whether 7 is traversed clockwise or anticlockwise). For the timelike curve to be 
closed {At = 0) we then see that 7 must be traversed clockwise (so that A < 0). 
Moreover, since we must necessarily have v < 1, the length / of 7 satisfies 

(6) 1<V2\A\. 

Given \A\, we known that / is minimized when 7 is a circle |Oss78[ IHHM99) . Now 
a circle of radius r in the hyperbolic plane satisfies 

I — 27rsinhr and \A\ — 27r(coshr — 1), 



■^Notice that none of adimensional quantities which we will compute - e.g. angles and velocities, 
including the total integrated acceleration - depend on the choice of units. 

^There is a sign difference corresponding to a different choice of the direction of rotation of the 
fluid elements, which we take to be anti-clockwise. 
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and so ([6]) is satisfied by a circle of radius r if 

sinhr < V2(coshr — 1), 

that is, if r is bigger than the positive root of 

sinh^ r = 2(coshr - 1)^ cosh^ r - 4 cosh r + 3 0, 

which is r = arccoshS. We conclude that ^ can only hold if |A| > 4tt. In 
other words, the projection 7 of a closed timelike curve must enclose a sufficiently 
large area. Notice that ii \A\ > An then there is indeed a closed timelike curve 
whose projection encloses this area - just choose 7 to be a circle and u < 1 an 
appropriate constant. Also, the isoperimetric property of the circle implies that the 
the projection of a closed timelike curve must have length / > 2tt sinh(arccosh(3)) = 
47rV2- Moreover, we see from ([5]) that the maximum velocity of the closed timelike 
curve with respect to the matter, Umax, must satisfy 

/- , Z ^ I sinhr 

Vmax V2\A\ %/2(coshr-l) 2 

again by the isoperimetric property of the circle (an equivalent observation was 
made in ;CGL83;). 

Choosing the natural orthonormal frame 

1^0 — — 



dt 

_ \__d_ I V2(coshr~l) d 

^ sinh r sinh r dt 



with dual coframe 



= d< + V2(l - cosh r)civ' 
cj'' = dr 
bj^ = smhr dip 

one readily obtains from the first Cartan structure equations 

dw" + w"^ A w'^ = 
the nonvanishing connection forms 

(7) 



The second Cartan structure equations 



= -ujf^ = - cotanhrw"^ 



then yield the nonvanishing curvature forms 



r 


_ or _ _ 

— Si Q — 


30; 


Aw'' 




— S2 Q — 


1, ,0 


A 




r 


1/ 
2^^ 





From the general expression 
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we obtain the independent components of the Riemann tensor in this orthonormal 
frame, 

(8) RorOr = RoipOip = Rriprip = 

Consider now a timeUke curve whose projection on the hyperbohc plane is a 
circle of radius r. Because of the space-time homogeneity of the Godel universe, we 
may assume that the circle is centered at the origin r = of our coordinate system. 
If the circle is traversed clockwise with constant velocity v = tanhw with respect 
to the matter then the four-velocity of the timelike curve is 

U = cosh u Eq — sinh u , 

and hence its covariant acceleration is 

(9) VuU = coshuuj%{U)Ea ~ sinh uw"^ 

= sinhwfvzcoshu — sinhucotanhr)^^^- 



In particular, the curve is a geodesic if and only if w = v^tanhr. The space-time 
homogeneity of the Godel universe implies that all timelike geodesies of the Godel 
universe are of this form. Notice that the radius of the circle satisfies 

tanhr < 4^ r < arccosh(\/2), 

and so no timelike geodesic can be closed. 

3. The Jacobi equation 

The minimum conditions ([T]) , ([2]) on a free-fall arc reduce to the Jacobi equation 
for a Jacobi field orthogonal to the four- velocity. If we choose coordinates (r, 1^9) 
centered at the center of the free-fall arc, we have 

U = cosh u Eq — sinh u E^ 

with tanhu = \/2 tanhr, and so 

P = F'Er + P'^E^, 

where 

E^ = — sinh u Eq + cosh u E^ 

is the unit vector field orthogonal to U and Er pointing in the same direction as 
Eip. Using the connection forms ([T]) and tanhu = \/2 tanhr, one readily checks 
that ^ 

{VuEr = —^E^ 
VljEip — -^Er 

for a geodesic, and so 

VuP = (^"^ 

whence 

Vc/Vt/P =(p'' + a/2P'^ - Er+iP"^ - V2P'' - l-P"" ] E^ 




2 \ 2 



The Jacobi equation, 
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has three components. The component along U, 

is trivially satisfied. The component along Er, 

P'' + V2P'^ ^ ^■^'^ " --Rroro cosh^ u P'' - i?,.^^^ sinh^ u P'' 
and the component along E, 



coshu ^P^ - %/2P'' - ip*^^ = -iJ^o^oCOsh^uP'^ + P^oo>pcoshu sinhuP", 
form the system of linear ODEs 

J P'^ + V2P'^ = - sinh^ u P'' 
\p^ - V2P'' = 
(where we used ([5])). This system is readily solved to yield 

(pr ^ 4_5cos(tJT) - Csin(a;r) 

\p'^ = Ail - + ^ sin(wT) + ^ cos(a;r) + £> 



where A,B,C,D e M are integrations constants and uj — \/ cosh^ m + 1. Notice 
that the period of the projection of the timelike geodesic on the hyperbolic plane 
is precisely 

27r sinh r 27r tanh r \/2-k tanh u 27r 

sinhw " sinhuVl-tanhV " sinhT/y/f - i tanh^ ~ ^ ' 

that is, Lo is the angular frequency with respect to the proper time r. 

4. Optimal circular accelerated arcs 

For an optimal circular accelerated arc the covariant acceleration Q is propor- 
tional to Er, and so P = ±Er- Using the connection forms ([7]) one readily checks 
that 

V2 . (V2 . \ 

VijEr — sinhwi?o + I coshw — cotanhrsinhu 1 E^, 

and so 

(V[/P;^)(V;7P^) = - — v^ sinh u cosh u cotanh r + cotanh^ rsinh^ u. 
From (HI) we have 

= s'mh^ u cosh^ u — 2-\/2 sinh u cosh u cotanh r + sinh^ u cotanh^ , 
and from ([8]) 

K ~ —Rai3f_iuP"U^ P'^U'^ = — cosh^ wPj-oro— sinh^ uRr,^r^ — ^2 ^"^^^^ " 2 ^^'^^^ 
Therefore we can write ^ as 



, 2 %/2 , 1 

tanh r tanhr H ^ = 0, 

u 1 + 
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where v ~ tanh u is the constant velocity of the optimal accelerated arc with respect 
to the matter. This can be solved with respect to tanhr to yield 



The ± sign determines whether VijU points in the same direction as Er or in the 
opposite direction. Indeed, from ([9]) we see that this happens if and only if 



a/2 — > <:=^ \/2vtanhr > <^ 1 ± \ ^ > <^ ± , > -1. 

tanh r ^ I + ^1 _ ^4 

It is easy to see that as v varies from to 1 the — branch of ([TU]) increases from 
to whereas the + branch decreases from +oo to Since we must have 
tanhr < 1, the + branch is only defined for w > iiq — 0.916126. 

Since tanh(arccosh(3)) = > we see that only the + branch of ([TOl) can 
contain an optimal circular accelerated arc which is a closed timelike curve. Such 
a curve must satisfy 

I- . f dl I—, , 27r sinhr 
\J2A + (p — ^ ^ -27rV2(coshr - 1) H sinhr ^O^v 



V V 



\/2(coshr- 1)' 



Therefore we must have 

1 — 2(coshr — 1)^ — sinh^ r cosh r — 3 



1 + u^ 2(coshr - 1)2 + sinh^ r 3coshr-l 

and so from (1101) 



= 1 + \ h: ^ r ^ 



sinh^ r / coshr — 3 1 / coshr — 3 



cosh r (cosh r — 1) V3coshr~l coshr V3coshr— 1 

Squaring this equation we obtain 

(11) cosh^ r — 3 cosh^ r — 3 coshr + 1 = 0, 

and, since x = — 1 is an obvious root of the polynomial x'^ — 32;^ — 3x + 1 , we readily 
obtain r = arccosh(2 + a/3) as the only real solution of ([TT|) . 

This optimal closed accelerated arc is precisely Malament's conjectured optimal 
trajectory. Its total integrated acceleration is 

27r sinh r ■ , , fK ^ . i i n 27r sinh r 

TA =^ a — — = smhw(v2coshM — smhucotanhr) — — , 

sinh u sinh u 

or, since 

tanhu = _ <^ coshusinhr = a/2 sinh u( coshr — 1), 

\/2(coshr-l) 

we have 

TA = 27r sinh w (coshr - 2). 
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Using 



'coshr = 2 + V3 



sinh r = \/ cosh'^ r — 1 = -y/ 6 + 4^3 



sinh r 



\/2(cosh r — 1^ 

cosh u ~ ^ 



2 



sinh u = 



= V4 + 273 
= V/3T2VI 



9 + 6\/3~ 27.6691, 



we obtain 

TA = 2tt 
as advertised in the introduction. 

The optimality of this curve can be partially checked as follows. Using the 
formulae above, it is not hard to show that 



TA = 27r(coshr- 2] 



cosh r + 1 



cosh r — 3 

for closed timelike curves with constant velocity with respect to the matter whose 
projection on the hyperbolic plane is a circle of radius r. It is then easy to show 
that Malament's curve does indeed minimize TA among these curves. 
Notice that the proper time necessary to traverse Malament's curve is 

27r sinh r r- 

—— = 2\/27r, 

smhu 

which in our units can be written as 

2V27r _ Fk 

~ V p' 

Therefore if the matter density of the Godel universe is comparable to the matter 
density of our universe then the proper time for this curve is comparable to the 
Hubble time for our universe, about 10^° years. The proper acceleration will be 
correspondingly low. 

5. The optimal trajectory 

We now present our candidate for the optimal trajectory and give evidence for its 
optimality. Our simplifying assumption is that this trajectory contains a circular 
accelerated arc of radius R> 0. Assuming, by analogy with Malament's trajectory, 
that the acceleration points outwards, we have from ([TU]) that the velocity V of this 
arc with respect to the matter satisfies 



(12) 



tanh R 



1 



V2V 



1 



1 - U2 

1 + 



Because we are simply looking for the closed timelike curve which minimizes the 
total integrated acceleration, we do not specify the initial and final four- velocities, 
and so, as stated in section [TJ P must vanish at the endpoints. Since p = 1 on 
accelerated arcs and instantaneous accelerations, we conclude that these endpoints 
must be on free-fall arcs. Moreover, we have P — Er throughout the accelerated 
arc, and so if there are instantaneous accelerations at its endpoints these must be in 
the outward radial direction. The simplest candidate curve of this kind is an initial 
free-fall arc, followed by the circular accelerated arc, followed by another free-fall 
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arc, possibly with instantaneous accelerations at the transitions. By simplicity we 
assume that the final free-fall arc mirrors the initial one, both having radius r > 
and velocity v — v^tanhr with respect to the matter. The projection of this 
candidate curve on the hyperbolic plane will then be an arc of a large circle of 
radius R, and two arcs of smaller circles of radius r, as shown in Figure [2] (where 
we use circles and lines to represent circles and geodesies of the hyperbolic plane) . 




Figure 2. Projection of the optimal trajectory on the hyperbolic 
plane (circles and lines represent circles and geodesies of the hy- 
perbolic plane). 

Each trajectory of this kind is determined by 3 parameters, which we will take 
to be the velocity V of the accelerated arc with respect to the matter, the boost 
parameter Am of the instantaneous accelerations, and the angle 9 in Figured This 
immediately fixes R through p^ . On the other hand, if we use coordinates centered 
at the center of the accelerated arc then the four- velocity of this arc is 

U-=TEo-VrE^, 

where 

If we apply an instantaneous acceleration in the radial direction with boost param- 
eter Au then the four- velocity jumps to 

U+ = cosh{Au)U^ + smh{Au)Er 

= rcosh(AM)£'o + smh{Au)Er - FF cosh(Au)£'^. 
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On the other hand, wc must also have 

C/+ = jEq + 117 sin aEj, — vj cos aE^p, 

where v is the velocity of the final (hence also the initial) free-fall arc with respect 
to the matter, 

1 

7 = 



and a is the angle between the two radii depicted in Figure[2] (notice that the radius 
of a circle in the hyperbolic plane is orthogonal to the circle, and so a is also the 
angle between the two tangents). Therefore 



ysinh^(Au) + 

7 = r cosh(Aw) ^ V 



cosh(AM) 
Dve 

tanh(Au) 



Notice that this fixes r = arctanh (^^^^ ■ Moreover, we have 

tan a= yj. ^ 

and so a is also fixed. Now if x is the distance between the center of the accelerated 
arc and the center of the free-fall arc, we have from the hyperbolic law of cosines 
|And07) 

cosh X — cosh r cosh R — sinh r sinh R cos a 

(see Figure [2]), which fixes x. From the hyperbolic law of sines we have 

sinh r 
sm p — ^— — sm a 
smh X 

and 

sinhi? 
sine = — — — sma, 
smh a; 

which fixes the angles /3 and s depicted in Figure [2] Similarly, if y is the distance 
between the center of the accelerated arc and the outermost intersection of the 
free-fall arcs, we have from the hyperbolic law of cosines 

cosh y = cosh r cosh x — sinh r sinh x cos 9 

(see Figure [5]), which fixes y. From the hyperbolic law of sines we have 

sinhr 
sm (fi — -— — sin Q 
smhy 

and 

sinh a; . 
sm?/' = -— — smb', 
smhy 

which fixes the final two angles Lp and -0 depicted in Figure [2] 

The three parameters V , Am and d are not independent, as we are looking for 
a closed timelike curve. Equation ([S]) implies that the area A enclosed by the 
projection curve in the hyperbolic plane must satisfy 

(13) A^^ + ^, 
where 

(14) L = 2(7r-/3-(^) sinhi? 
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is the length of the accelerated arc, and 

(15) ; = 2(27r-£-6i)sinhr 

is the sum of the lengths of the free-fall arcs. Using the Gauss-Bonnet formula for 
the area of a hyperbolic triangle jAndOTj , we can write 

(16) A^2{n- 13- ip){coshR - 1) + 2(27r - e - 6l)(coshr - 1) 

+ 2{Tr - a - P - e) + 2{Tr - ip - e - Tp). 

From ([9]), the acceleration the accelerated arc is 

a^Vr^{V2~V cotanh R) , 
and the proper time necessary to traverse it is 

2{Tr- P-ip) srnhR 

vf ' 

Therefore the trajectory's total integrated acceleration (including the two instan- 
taneous accelerations) is 

(17) TA = 2r{n ~ P- v3)(\/2sinhi? - V^coshi?^) + 2Au. 

Notice that the angle between the projections of the two free-fall arcs is 21/). This 
means that if {Fq, Fr, F^} is the orthonormal frame associated to coordinates (r, ip) 
centered at the center of the final free-fall arc, so that the final four- velocity is given 

by 

Ui = 7F0 - vjF^, 

then the initial four- velocity is 

Uq — jFq + vj sm(2'ip)Fr — vj cos{2ip)F^. 

Therefore the hyperbolic angle Am' between the corresponding four-velocities sat- 
isfies 

-coshAu' = (C/o,C/i) = -7^ + w^7^cos(2i/'), 

where (•, •) is the metric, and so the integrated acceleration necessary to boost the 
final four- velocity to the initial four- velocity is 

Au' — arccosh [ — j \ ^ 

\ \ — J 

If one wishes to regard the candidate curve as periodic then one must add this value 
of Au' to p7)) . 

We are now ready to search for the optimal curve. Because of the complexity of 
the formulae above, a numerical search seems to be inevitable. We use a shooting 
algorithm as follows: starting with arbitrary values for Am, 9 and V ^ we determine 
the remaining parameters R,r,x,y,a, I3,e,(p,tp. We then use p^ . and 
to compute L, I and A, and check whether holds; if not, we modify the value 
of V until it does. At this point we have a closed timelike curve and can use pT)) 
to compute its total integrated acceleration. Using this algorithm, we vary Au 
(keeping 9 fixed) until TA is minimized, and then vary 9 (keeping Au fixed) until 
TA is again minimized. Repeating these iterations the method quickly converges to 
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(19) 



yield a closed timelike curve whose total integrated acceleration is minimum among 
all curves of this type. The solution is determined by 

{V ~ 0.930597 
Au~ 1.16765 
e ~ 1.64430 

We notice in passing that the velocity of the accelerated arc is very close to, but 
slightly smaller than, the velocity of Malament's curve, which is 

~ 0.930605. 

The remaining parameters of the optimal trajectory are 

'v ~ 0.978202 
r ~ 0.851195 
R ~ 1.99194 
a ~ 0.313261 
X ~ 1.25222 
P ~ 0.184815 
e ~ 2.37999 
y ~ 1.66316 
V? ~ 0.385023 
~ 0.681328 

yielding the total integrated acceleration 

TA ~ 24.9927. 

If we boost the final four-velocity to the initial four- velocity to make the trajectory 
periodic then the total integrated acceleration becomes 

TA ~ 28.6085, 

in agreement with Malament's conjecture. 

We now present evidence that this is indeed the optimal solution we seek. To 
avoid confusion, we keep {Eq, Er, E^p} as the orthonormal frame associated to co- 
ordinates {r,Lp) centered at the center of the accelerated arc, and let {Fg, ^V, Fi^} 
be the orthonormal frame associated to coordinates (r, tp) centered at the center of 
the final free-fall arc (notice that i?o = ^o). Recall from section [3] that on the final 
free-fall arc the primer has components 

-^ + B cos(a;r) - C sin(wT) 
Ail-^)T + ^ sin(a;r) + ^ cos(wt) + D 



(20) 



in the orthonormal frame {[/, -Fi-, F^}, where 

U = 7F0 - v-fF^, F^ = -v-fFo + jF^, 
A, B,C, D G M. are integration constants and 
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If for simplicity we set r = as the instant of transition from the accelerated arc 
to the free-fall arc, we have 



(21) 



y pip _ V2C _|_ £) 

at that instant. Moreover, it is easy to check that 

f f"- = -Cui 



(22) 
and 
(23) 



at the same instant. 

Now, as we saw above, the instantaneous acceleration boosts the four-velocity 
from 

U- = TEo - VTE^, 

to 

U+ = r cosh(A«)i;o + smh{Au)Er - yrcosh(Au)i;^. 

Consequently the new tangential direction is given by (minus) the normalized spa- 
tial velocity, 



y^sinh^(Au) + V^T^ cosh^(Au) ' 

and consequently the new radial direction is given by 

VT cosh(Au)£r + smh{Au)Ey, 
tr = I 

Y'sinh^(Au) + V^V^ cosh^(Au) 

On the other hand, the primer is boosted from 

P_ = Er 

to 

P+ = sinh(Att)C/_ + cosh{/\u)Er 

= r sinh( Au)£;o + cosh(AM)_E;r - VT s\n\i{ Au)E^. 
Therefore the components of the primer at the beginning of the free-fall arc are 

(24) P^ = (P+,F,)= 



/sinh2(Aw) + V^T^ cosh2(Au) 
and 

. , ~,„ / ~ \ 1 , , — rsinh(A'u) 

(25) P^ = (P+,f^) = -(P+,f^)= ^ ' , 

^ ' ^ ^sinh2(Au)-hy2r2cosh2(Au) 

where we used 

= v^U+ + jF^ ^F^ = -F^ - vU+ 

and 7 = rcosh(Au). 
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Recall that we must have P = at the endpoint of the final free-fall arc. Now 
from Figure [2] it is clear that lot = 27r — e — at this endpoint. From (|20|). ([2l]) . 
and ^ we then have 

(-:4a + B = - 

\/2(^ _!_£) — -rsinh(Ati) 

-^yl + cos(27r -e-e)B - sin(27r - e - 6l)C = 

. (1 - ^''"J"^ ^ + ^ sin(27r -e-9)B + ^ cos(27r - e - 6l)C + D = 

This can be regarded as a linear system in the integration constants A, B,C, D ^ R 
which determine the primer in the final free-fall arc. Using the numerical data in 
and (HH) we obtain 

A ~ 1.04251 
B ~ 0.600629 
C ~ -0.572707 
D ~ -0.677195 

Now on the accelerated arc p is identically 1. Since, as stated in section [1] p is a 
function of class C^, we must have p = p = at the instantaneous acceleration. 
Given that 

we have 

pp P'-p'- + P'fip'fi 

and 

2*2 

p'^ + PP^ (/"■) + (p'^) + P^'P'' + pfp'f. 

Using (pTjl , , and the numerical results in I^E^ we then obtain 

-8.71853 X 10"*^ 

and 

p2 -I- ~ 2.79359 X 10-^ 

at the instantaneous acceleration, which is very close to zero indeed. Notice that 
although p = 1 at the instantaneous acceleration by construction, the values of p 
and p on the final free-fall arc arise from the Jacobi equation, and are in principle 
completely independent. The fact that they both essentially vanish at the instanta- 
neous acceleration is powerful evidence in favor of the optimality of our candidate 
curve. 

6. Conclusion 

In this paper we presented a candidate for the closed timelike curve in the Godel 
universe which minimizes the total integrated acceleration. The projection of the 
candidate curve on the hyperbolic plane representing the spatial part of the Godel 
universe is depicted in Figure [2l and is traversed clockwise. It starts with a free-fall 
arc of velocity v ~ 0.978202 with respect to the matter, whose projection on the 
hyperbolic plane a 129° arc of a circle of radius r ~ 0.851195 in our units. This 
free-fall arc ends with an instantaneous acceleration of total integrated acceleration 
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Au ~ 1.16765, which, as seen by the matter, rotates the spatial velocity anti- 
clockwise by 18° and reduces its magnitude to V ~ 0.930597. The candidate curve 
then coincides with a circular accelerated arc of constant velocity with respect to 
the matter, whose projection on the hyperbolic plane is a 295° arc of a circle radius 
R ~ 1.99194. This sequence of events is reversed at the end of the accelerated 
arc: there is a new instantaneous acceleration of total integrated acceleration Au ~ 
1.16765, which, as seen by the matter, rotates the spatial velocity anti-clockwise by 
18° and incrciases its magnitude back to w ~ 0.978202, after which tlic;re is a new 
free-fall arc whose projection on the hyperbolic plane is another 129° arc of a circle 
of radius r ~ 0.851195. 

Notice that while our candidate curve does satisfy the minimum conditions, it 
is possible that it is just a local minimum, or even a saddle point. The question 
of whether it is indeed the minimum total integrated acceleration closed timelike 
curve in the Godel universe will have to be settled by different methods. 
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